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7—1 -Abstract 

o 

CS| 'Since the Voyager fly-bys of Uranus and Neptune, improved gravity field data have been derived from long-term obser- 
vations of the planets' satellite motions, and modified shape and solid-body rotation periods were suggested. A faster 
rotation period (— 40min) for Uranus and a slower rotation period (+lh20) of Neptune compared to the Voyager data 
were found to minimize the dynamical heights and wind speeds. We apply the improved gravity data, the modified 
shape and rotation data, and the physical LM-R equation of state to compute adiabatic three-layer structure models, 
.where rocks are confined to the core, and homogeneous thermal evolution models of Uranus and Neptune. We present 
r— I'the full range of structure models for both the Voyager and the modified shape and rotation data. In contrast to pre- 
.vious studies based solely on the Voyager data or on empirical EOS, we find that Uranus and Neptune may differ to 
M-l 'an observationally significant level in their atmospheric heavy element mass fraction Z\ and nondimensional moment 
^ !of inertia, A. For Uranus, we find Z\ < 8% and A = 0.2224(1), while for Neptune Z\ < 65% and A = 0.2555(2) when 
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O-i applying the modified shape and rotation data, while for the unmodified data we compute Z\ < 17% and A = 0.230(1) 
for Uranus and Z\ < 54% and A = 0.2410(8) for Neptune. In each of these cases, solar metallicity models [Z\ = 0.015) 
are still possible. The cooling times obtained for each planet are similar to recent calculations with the Voyager rotation 
periods: Neptune's luminosity can be explained by assuming an adiabatic interior while Uranus cools far too slowly. 
More accurate determinations of these planets' gravity fields, shapes, rotation periods, atmospheric heavy element abun- 
dances, and intrinsic luminosities are essential for improving our understanding of the internal structure and evolution 
of icy planets. 
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1. Introduction 

The outer planets Uranus and Neptune are mysterious in 
many ways. While their names 'ice giants ' suggest a com- 
position of predominantly volatiles in ice phases such as 
water, methane, and ammonia ice, interior models instead 
predict a warm interior devoid of solid ices (Chau et al., 
2011: Redmer et al., 2011). Structure models of Uranus 
and Neptune are generally in agreement in predicting a 
small rock core, a deep interior of more than 70% heavy 
elements, and a significantly less enriched outer envelope 
with a transition at about 70% of the radius for both plan- 
ets (Hubbard et al., 1995; Fortncy and Ncttclmann, 2010; 
Helled et al., 2011). However, they also agree in failing 
to explain Uranus' measured intrinsic luminosity (Podolak 
et al., 1991). Uranus' low luminosity is a riddle; even more 
so as the corresponding models for Neptune nowadays can 
reproduce its measured luminosity (Fortney et al., 2011) 
indicating that the ice giants are not so similar to each 
other as previously thought (Podolak et al., 1995). 

Contrary to the efforts that have been made to mea- 



sure the atmospheric composition, the Voyager 2 radio oc- 
cupation data and ground-based observational microwave 
data have actually raised more questions than they were 
intended to solve. In both planetary atmospheres for in- 
stance, D/H appears enriched over the protosolar values 
(de Bergh et al., 1986, 1990) suggesting particle trans- 
port between the atmosphere and an ice-rich deep in- 
terior (Hubbard et al., 1995), while the inferred helium 
abundance is consistent with the protosolar value and 
thus would suggest the opposite (i.e., inefficient particle 
transport) if hydrogen and helium occur in the deep inte- 
rior where they might undergo phase separation (Hubbard 
et al., 1995; Podolak et al., 1995). Moreover, for the inte- 
rior of Uranus we cannot rule out a primary composition of 
silicates dissolved into hydrogen-helium envelopes as that 
might sufficiently accelerate Uranus' cooling time (Hub- 
bard and MacFarlane, 1980). We just do not know how 
similar Uranus and Neptune really are in terms of compo- 
sition and structure (Podolak et al., 2000). 

In addition, Uranus and Neptune have complex multi- 
polar magnetic fields, and appear to have stronger atmo- 
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spheric winds than Jupiter and Saturn, if the Voyager ro- 
tation periods are applied. Besides the natural desire to 
understand the giant planets in our solar system, learning 
about the 'ice giants' is important for the classification of 
extrasolar planets with similar masses and sizes, as many 
ones are observed (Borucki et al., 2011), although not yet 
on similar orbital distances (Kane, 2011). 

Information on Uranus and Neptune interiors are typi- 
cally derived from theoretical models which are designed 
to fit the observed physical data of the planets, such as 
their gravitational fields, masses, internal rotation, and 
radii. The physical data available for Uranus and Nep- 
tune are rather limited. In particular, the low-order grav- 
itational harmonics J2 and J4 of Neptune have significant 
error bars, Neptune's equatorial 1-bar radius is actually 
not measured, and the shapes (flattening) of both Uranus 
and Neptune are not well known. Voyager 2 provided only 
one occultation radius at the 1 bar pressure level for each 
planet (see Helled et al., 2010, for details). Although stel- 
lar occultations do provide information on the planetary 
oblateness, the shape is inferred for upper atmosphere (mi- 
crobar pressure levels) and it is unclear whether this shape 
is consistent with the shape at the 1 bar pressure level. 
Helled et al. (2010) have shown that minimization of wind 
velocities or dynamic heights of the 1 bar isosurfaces (An- 
derson and Schubert, 2007; Helled et al., 2009) constrained 
by Voyager 2 occultation radii and gravitational coeffi- 
cients of the planets, leads to modified solid-body rotation 
periods of 16h 34m for Uranus and 17h 27m for Neptune, 
which is 40 min shorter than the Voyager rotation period 
for Uranus, and lh20 longer for Neptune. Helled et al. 
(2010) also state that both planets may have different ro- 
tation periods than the ones derived by the minimization 
method, and in addition, could be rotating differentially 
on cylinders. Non-solid body rotation can lead to a change 
in the calculated gravitational moments (Hubbard et al., 
1991). Our results are valid as long as differential rota- 
tion in Uranus and Neptune is "shallow", i.e., the region 
of differential rotation consists of a negligible mass, which 
is a preferred solution for Uranus and Neptune (Y. Kaspi, 
priv. comm.). However, the suggested solid-body rotation 
periods better match the measured radii of the planets 
and result in more moderate atmosphere dynamics with 
wind velocities of ~ 150 m s _1 for both planets. Based on 
these suggested rotation rates and occultation radii, mod- 
ified shapes of the planets were derived. Since the shapes 
of Uranus and Neptune are not well constrained, and in 
addition, their rotation profiles must not be that of a solid- 
body, the modified solid-body rotation periods and shapes 
can be used to demonstrate the sensitivity of the interior 
models to the uncertainties in these physical properties. 

To obtain pressure-density relations for the interior of 
Uranus and Neptune, three different methods have been 
invoked. Podolak et al. (1995) apply physical equations 
of state of hydrogen, helium, the ices H2O, H2S, CH4, 
and NH3, and rocks, assuming sharp transitions between 
a gaseous outer envelope, an inner icy envelope, and a 



rock core. As this assumption by itself restricts the possi- 
ble internal density distributions, Marley et al. (1995) and 
Podolak et al. (2000) created random density distributions 
which, when reproducing the observed gravity field, could 
then be evaluated according to pressure-density relations 
of likely materials. Another possibility is to use an an- 
alytic function with sufficiently many free coefficients to 
adjust to the given constraints (Helled et al., 2010). The 
resulting pressure-density relation can be considered an 
empirical EOS. The advantage of the two latter methods 
is that they can allow for, and one can constrain the loca- 
tions of, continuous (rather than sharp) density gradients. 
On the other hand, the resulting random or fitted density 
distributions can be unrealistic in a sense of not represent- 
ing any composition of real materials. We here use the 
physical EOS LM-REOS for H, He, and H 2 (Nettelmann 
et al., 2008) for the envelope material and the rock EOS by 
Hubbard and Marley (1989) for the core. This combina- 
tion was also applied to the Uranus and Neptune models 
in Fortney and Nettelmann (2010) (hereafter FN10). 

In this paper, we use the suggested modified periods and 
equatorial radii and the physical LM-R EOS to compute 
new three-layer interior and evolution models of Uranus 
and Neptune. In § 2 we describe our method of model 
construction and the observational data used. In § 3 we 
compare the resulting models obtained with the modified 
data to those with the Voyager data. The results are dis- 
cussed and summarized in § 4. 

2. Method of Structure and Evolution Modeling 

Our Uranus and Neptune models consist of three lay- 
ers; a rocky core surrounded by two adiabatic and homo- 
geneous envelopes of hydrogen, helium, and water. The 
interior model uses the following physical data: the plane- 
tary total mass M p , the equatorial radius R eq , the temper- 
ature at the 1 bar pressure level T\ which determines the 
internal adiabat, the rotation rate ui = 2tt/P with P being 
the solid-body rotation period, and the measured gravita- 
tional coefficients J2, J4. The physical parameters used 
in the models are listed in Table 1. For simplicity, it is 
assumed that both Uranus and Neptune arc adiabatic and 
convective, although Uranus may not have an adiabatic 
temperature gradient (Podolak et al., 1991). The resulting 
model parameters of interest here are the mass fraction of 
heavy elements in the inner envelope, Z2, and in the outer 
envelope, Z\, where we restrict Z\ to > 1 Zq, Zq = 0.015 
being the solar metallicity (Lodders, 2003), the core mass, 
M c , and the possible transition pressure Pi_2 between the 
envelopes. The resulting values Z\, Z2, and M c are ob- 
tained by an iterative scheme that fits J2, J4 and the con- 
dition of mass conservation m = at the center r = 0. A 
constant Hc:H mass ratio of 0.275:0.725 is assumed in the 
envelopes. The method used to calculate interior models 
is the same as in (Nettelmann et al. 2008, FN10) and the 
thermal evolution procedure the same as in FN10. The 
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progress over that work is mainly the application of im- 
proved observational data. Recent work by Jacobson and 
collaborators provided new estimates for the gravitational 
harmonics of both Uranus and Neptune using Earth-based 
astrometry and observations acquired with the Voyager 
spacecraft. Their more accurate determinations included 
more astrometric observations over the years, a change in 
reference frame to the International Celestial Reference 
Frame, a modification in the numerically integrated Pro- 
teus orbit (for Neptune), and an improved data processing 
procedure (see Jacobson (2009) for details). As in Helled 
et al. (2009, 2010) we here use the Jacobson (2007, 2009) 
data 1 for the gravitational moments J 2 and J4. These are 
provided at reference equatorial radii as given in Table 1. 
The gravitational moments at the equatorial 1-bar radius 
to be fitted by our models are then obtained from the mea- 
sured values (primed quantities) by scaling according to 
the relation P^™ J2„ = R'cqJ^n- Figure 1 shows these im- 
proved data in comparison with the former Voyager grav- 
ity data, and the relative observational uncertainties for 
all four outer planets. 

[Table 1] 
[Figure 1] 

For the thermal evolution we make the standard as- 
sumption that the planet cools down over time by the 
release of gravitational and internal energy remnant from 
their formation, where the heat from the interior is trans- 
ported by convection along an adiabatic temperature gra- 
dient. Thus our assumption of different heavy element 
mass fractions in the two envelopes, (Z\ ^ Z2) implies 
that one compositional gradient exists that inhibits con- 
vection across the layer boundary, but not the heat flow. 

With the cooling time r we denote the time the planet 
needs to cool down from an arbitrarily hot initial state to 
its current effective temperature T c g. To calculate r, we 
pick a certain structure model (meaning a set of stationary 
parameters {Z±, Z2, mi-2, M c }, where mi_2 = m(Pi_2)), 
and calculate ~ 60 internal profiles with stepwise increased 
T\ values, that represent the planet at earlier times. We 
then integrate the cooling equation L e s — L cq = —dE^/dt 
over time, where L e ff is the measured infrared luminos- 
ity, -L oq the insolation, and dEi nt the lost intrinsic energy 
during the time interval dt. To close the relations, we 
use the model atmosphere approximation by Guillot et al. 
(1995), and the Stefan-Boltzmann law to relate T\ to L e ff- 
While the underlying grid for non-irradiated, solar metal- 
licity model atmospheres by Graboske et al. (1975) may 
appear inappropriate for application to Uranus and Nep- 
tune despite its former application (Hubbard and Mac- 
Farlane, 1980, FN10), the induced uncertainty in the re- 
sulting cooling times of Uranus and Neptune compared to 
those based on a more self-consistent grid for Uranus and 
Neptune has been found to be ~ 0.2 Gyr only (Fortney 
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et al., 2011). Both grids neglect the presence of clouds 
and hazes which are thought to be present in the tropo- 
spheres (Gautier et al., 1995). We note that clouds are 
just at the beginning of being incorporated into climate 
models for the Earth (Dessler, 2010). Therefore, the full 
uncertainty on the cooling time from the application of 
the model atmosphere grid is not known but future work 
on cloudy skies is not expected to significantly alter the 
main findings of this work as we are searching for differ- 
ences between Uranus and Neptune. This also justifies the 
ncglection of the time-dependence of the insolation, of the 
rotation rate due to angular momentum conservation, and 
of the energy of rotation. Their effects on the cooling time 
have been shown to nearly compensate each other in case 
of Jupiter (Nettelmann et al., 2012). 

3. Results 

3.1. Models with the Voyager shape and rotation data 

Figure 2 shows our results for the envelope metallicitics 
of three-layer Uranus and Neptune models with the Voy- 
ager shape and rotation data. There are no solutions with 
Z\ > Zq outside the respective boxes, as either the core 
mass becomes zero, or | J4I becomes too low or too high. 

Models with the same transition pressure show a lin- 
ear dependence between Z\ and Z2, because to first order, 
the resulting decrease in the outer envelope's mass {dM\) 
is proportional to the imposed reduction —dZ\. Mass 
conservation then requires dM2 + dM c = —dMi, where 
dM2 ~ dZ2- Consequently, dZ2 ~ —dZ\. If mass is shifted 
downward {dM\ < 0), \J±\ must decrease. Thus a lower 
limit of Z\ at given P1-2 arises from the condition that 
I J4I must not become too small. Since M c decreases with 
increasing Z2 in standard three-layer models as known for 
Jupiter (Nettelmann et al., 2012), at very high Z2 values 
(> 0.95) the core mass decreases to zero before the lower 
observational limit of | J4I is reached. Thus the Neptune 
models along the upper right boundary in Fig. 2 are in fact 
two-layer models with no core. The internal adiabat then 
reaches temperatures up to 6600 K and pressures up to 
8 Mbar, while the lowest core-mantle boundary tempera- 
ture, T c (pressure, P c ) is 5200 K (5.5 Mbar), found for the 
model with the biggest rock core of mass M c = 3.7 M^. 

The described behavior is the same as found previously 
for Uranus and Neptune (see figure 3 in FN10). How- 
ever, due to the large observational error bars of J2 and 
J4 (Fig. 1) allowed for in FN10, the Neptune models were 
arbitrarily forced to have Z\ < 0.40 and Pi_2 > 0.15 Mbar 
in accordance with the results for Uranus. It is the much 
-reduced error bars of J2 and J4 that here allow us to 
present the full range of models. The biggest difference 
to the set of Neptune models in FN10 is the existence of 
models with pure water envelopes. Those occur for high 
outer envelope metallicitics Z\ > 0.5. Then, Pi_2 needs 
to be put deep inside the planet (P ~ 1 Mbar) to ensure 
that I J4I does not get too large. 
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Few Uranus models can be found for the improved er- 
ror bars of J2 and J4 (Fig. 1). They have a well de- 
fined Z-2 = 0.9 ± 0.02, a tightly constraint position of 
the layer boundary at 0.12 ± 0.04 Mbar corresponding to 
r = 0.79 ± 0.02 R v and m = 0.93 ± 0.02 M v , a rather low 
outer envelope metallicity 1 Zq < Z\ < 12 Zq, a rock core 
of at most 1 M®, and central temperatures T c ~ 6000 K. 
An example is model Ul (Table 2). In Z\-Zi space, the 
Uranus models are within the range of the Neptune mod- 
els. 

[Figure 2, top left] 
[Figure 3, top right] 

3.2. Models with the modified shape and rotation data 

The explanations to Fig. 2 in § 3.1 apply also to Fig. 3. 
Wc here describe the differences that arise from using the 
modified shape and rotation data. 

Most obvious, the sets of solutions shift into opposite 
directions: the Uranus models move to the lower right 
corner of Fig. 3 where Z2 gets higher and Z\ smaller, 
increasing the density difference between the envelopes. 
In contrast, the Neptune models stretch into the upper 
left corner where Z\ is above 60% and the density dif- 
ference to the deep envelope less pronounced. As a re- 
sult, the sets of solutions become disjunct. The Uranus 
solutions are tightly constrained to have a low outer enve- 
lope metallicity (< 0.1), a high inner envelope metallicity 
(> 0.9), a transition far out (at > 0.9 Mtj), and a small 
core (< 1.M®). An example is model U2, see Table 2. 
Similar Uranus and Neptune solutions are still possible, 
for instance with Z\ ~ 0.1 and Z^ ~ 0.88 in both plan- 
ets. However, this would require an envelope transition in 
Neptune at ~ 0.05 Mbar, relatively far out at 0.88 i?N in 
the neutral, molecular fluid part, which is not a preferred 
solution for Neptune. The new finding here is that two 
kinds of qualitatively different Neptune models are pos- 
sible that both differ in the atmospheric heavy element 
abundance to an obscrvationally significant level from the 
Uranus models. The first kind of Neptune models are char- 
acterized by a rather large core (M c ~ 3Mq) and a modest 
heavy element difference between the envelopes (changing 
from from 0.6 to 0.8). Alternatively, small core models are 
possible with pure water envelopes and a transition deep 
inside at ~ 0.6 i? N and ~ 0.6 M N . Models N2a and N2b 
are respective examples (Table 2). 

These changes can easily be explained by the new ro- 
tation rates. A slower rotation (Neptune) means a lower 
centrifugal force. Matter is then less strongly pushed to 
the outer region. If the J2, J4 to be fitted remain about 
the same, a higher metallicity in the outer part of the 
planet is required. Along the adiabat, a higher metal- 
licity leads to lower internal temperatures. Therefore, 
Neptune models with high outer envelope metallicity can 
become rather cold, with core-mantle boundary tempera- 
tures around 5000 K only. 



[Table 2] 

Example density profiles. The input and resulting param- 
eters of representative models arc given in Table 2 and 
their density and mass profiles shown in Fig. 4. The Uranus 
models clearly stand out by their big jump in density at 
~ 75% of the planet's radius, where the density falls down 
to 30% of the inner envelope boundary value. With the 
end of the inner envelope, 90% of the total mass is already 
accumulated. For Neptune, a model can be found (N2b) 
where the density difference is only 20% and thus more 
similar to the transition in standard three-layer Jupiter 
models (10%, Nettelmann et al., 2012) than to Uranus. For 
Neptune model N2a, the transition occurs deep inside at 
60% of the radius, where 57% of the total mass is accumu- 
lated, which would be more similar to standard three-layer 
models of Saturn (Nettelmann, 2009) than to any other 
outer planet. For both planets, a typical mean density in 
the inner envelope is 3 gcm~ 3 . This (arbitrary) density 
threshold occurs at about 0.5 R p (Fig. 4, left panel) and 
0.4 M p (Fig. 4, right panel), and corresponds to pressures 
of 2-3 Mbars and temperatures of 4000-5000 K (Fig. 5). 

[Figure 4] 
[Figure 5] 

Machine-readable tabulated interior profiles of mass, 
pressure, radius, temperature, and density of the models 
Ul, U2, Nl, and N2b are provided as supplemental online- 
material to this article. Table 3 show an example of such 
a table. 

[Table 3] 

Moment of Inertia. The above described dichotomy in 
the internal structures, which mainly refers to the outer 
50% of the planet's radius, maps onto different nondimen- 
sional moments of inertia 2 , A. For Uranus and Neptune, 
respectively, we calculate values of A = 0.2224(1) and 
A = 0.2555(2) when using the modified shape and rota- 
tion data, while in the unmodified case the respective val- 
ues arc A = 0.230(1) and A = 0.2410(8), where the number 
in parenthesis gives the uncertainty in the last digit, see 
also Figs. 2 and 3. 

3.3. Thermal evolution 

The cooling times for the representative models are 
given in Table 2. The modified rotation and shapes have 
essentially no effect on the cooling time. This is to be 
expected as the cooling time mainly depends on the spe- 
cific heat c v of the bulk material, while the composition 
(and hence c v ) is little influenced by the shape and rota- 
tion data. Interestingly, Neptune appears to be the outer 



2 A = // (M p iJj^ean ) , where / is the dimensional axial moment of 
inertia. 
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planet with the cooling time that best matches the age of 
the solar system, r Q = 4.56 Gyr. Under the assumption of 
adiabatic, homogeneous cooling, Saturn's theoretical cool- 
ing time is systematically too short by 2-2.5 Gyrs (Guillot 
et al., 1995; Fortney and Hubbard, 2003), while Jupiter's 
is often found to be a little too long (Saumon et al., 1992; 
Guillot et al., 1995; Fortney et al., 2011). We here re- 
cover the well-known far too long cooling time of Uranus 
(Fortney et al., 2011), indicating that some part of the in- 
terior is not in a state of efficient energy transport through 
vigorous convection (Podolak et al., 1991; Hubbard et al., 
1995). To estimate the size of the stable internal region 
in Uranus, we re-calculate the cooling time for the same 
quasi-adiabatic three-layer models as before but exclude a 
central mass, the cut-off mass m cut , from efficient cooling. 
The efficiency factor / quantifies the heat escape from the 
stable interior beneath m cut in comparison with the con- 
vective heat flux, so that / = implies zero heat flux and 
/ = 1, as well as m cut = 0, is the uninhibited case. In Fig. 6 
we present the cooling time of Uranus in dependence on 
m cut and /. 

[Figure 6] 

The larger m cu t , i.e. the smaller the convective outer re- 
gion, the shorter the cooling time. For a minimum cut-off 
mass of 0.45 My, t = r Q can be achieved. The larger the 
heat flux from the stable interior, the longer the cooling 
time, and thus the larger the required cut-off mass. A 
maximum value for / is obtained when m cut approaches 
M\j to realize r = r Q . While this maximum is / = 0.5, 
we consider small values / < 0.1 more realistic as the 
heat transport in a stable interior is predicted to be a tiny 
fraction of the convective heat flux. For Uranus, a lower 
limit for / can be the ratio F con d, a d/-Pconv ~ 0.02 , where 
^cond.ad is the conductive heat flux along an adiabatic gra- 
dient, and F conv the heat flux if Uranus would have cooled 
down to its present luminosity by large-scale convection 
(Podolak et al., 1995). Thus we conclude from figure 6 an 
ending of the stable interior in Uranus at 0.45—0.5 Mtj. 
For comparison, we have applied the same procedure to 
Neptune (model N2a). Neptune's cooling time appears 
weakly sensitive to the innermost 20% of Neptune's mass. 

4. Discussion and Conclusions 

4-.1. Bulk composition 

We had to make several simplifications to be able to cal- 
culate Uranus and Neptune models with state-of-the art 
methods. One such simplification is the representation of 
heavy elements in the envelopes by water and the confine- 
ment of rocks to the core. In real Uranus and Neptune, 
silicates may also occur in the envelopes. It is clear that 
our assumption of having no rocks in the envelopes leads 
to an overestimation of the envelope metallicities and the 
resulting ice to rock ratio (I:R). In addition, the smaller 



the core mass, the larger the I:R ratio. Models with I:R 
3> 2.7, the solar system value, may potentially invalidate 
this simplification and point to the presence of silicates in 
the envelopes. Example models are the Neptune models 
with pure water envelopes (N2b) and all our Uranus mod- 
els. Interestingly, some Neptune models, e.g. Nl and N2a, 
have a rather large core (~ 3M ffi ) with a reasonable overall 
I:R ratio of 1.5 times the solar value. On the other hand, 
our representation of ices by a pure water EOS likely over- 
stimates the density of the true mixture of ices, and thus 
somewhat underestimates the mass fraction of ices that 
would be composed of a mixture of H2O, CH4, and NH3. 
Unfortunately, high-quality EOS of light ices at pressures 
higher than 2 Mbar for planetary modeling are not yet 
available. 

Under the assumptions and simplifications of this work, 
the bulk mass of heavy elements is 12.5 M ffi for Uranus 
and 14-14.5 M® for the selected Neptune models. This is 
in good agreement with the empirical EOS based models 
by Helled et al. (2011), who can explain the polynomial 
density distributions of their Uranus (Neptune) models by 
~ 11-13 M ffi (~ 13-15 M & ) of heavy elements. 

4-2. Implications from a stable deep interior 

If some part of the interior, for Uranus possibly 0.45- 
0.5 My (§ 3.3), is stable to convection so that heat cannot 
escape efficiently from the region below, then the super- 
adiabaticity of the temperature gradient there can be non- 
negligible (Leconte and Chabrier, 2012). A warmer deep 
interior will require a lower particle number density in or- 
der to conserve the pressure gradient. Otherwise, the in- 
duced higher warm-dense-matter pressure would cause a 
larger planet radius. For a given composition in the in- 
ner envelope (e.g. the Z2 value of the adiabatic case), one 
might think of a lower particle number density to imply a 
lower envelope mean density (compared to the adiabatic 
case), resulting into a larger rock core mass to ensure mass 
conservation. However, the mean density in the inner en- 
velope is roughly constrained by the measured J2 value, 
see Fig. 4. Therefore, the mctallicity in the deep interior 
cannot have the same Zi value as in the adiabatic case. 
Indeed, if the stable region is caused by a compositional 
gradient (Hubbard et al., 1995), the deep interior would 
have an average mctallicity Z3 > Z2, which may lead to 
a smaller rock core. Therefore, predictions on a change of 
the rock core mass are not possible without more sophisti- 
cated models that take into account both a compositional 
gradient and super-adiabaticity. At the current level of 
our models, the net effect of a stable interior would be the 
introduction of a third, deep envelope with Z3 > Z2. In 
case of a significant super-adiabaticity in Uranus, this Z3 
might raise to 100%. Moreover, a H/Hc-frec deep inte- 
rior might then be possible even with a material density 
larger than that of water, i.e. with an ice-rock mixture. 
Thus, a strongly superadiabatic deep interior (below about 
0.5 My) may qualitatively allow for I:R ratios closer to the 
solar and at the same time for a H/He-free deep interior, 
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resembling a traditional protoplanetary core of a few M s 
that can accrete its gaseous envelope within a few million 
years (Hori and Ikoma, 2010). 

A stable deep (i.e., far below the outer/inner envelope 
boundary) interior offers an explanation for the measured 
atmospheric helium and deuterium abundances. Given 
that D/H can be enriched in the cold ices of the proto- 
solar nebula, the observed enhanced atmospheric D:H ra- 
tio could result from upward-transport of deuterium from 
an ice-rich interior (Gautier et al., 1995) that is located 
between the inner/outer envelope boundary and the sta- 
ble deep interior, where the magnetic field is believed to 
be generated, see § 4.3. Assuming that hydrogen and he- 
lium exist in the deep interior at Mbar pressures where 
hydrogen is metallic, Hubbard et al. (1995) suggest that 
helium could phase separate and rain down leading to a 
helium-depleted atmosphere. However, if the particle ex- 
change with the upper regions is suppressed in the deep 
interior, the rained out helium will not be replenished from 
the reservoir above so that the atmospheric Hc:H ratio re- 
mains protosolar as observed. However, we caution against 
taking the observed D and He abundances for clear indica- 
tions of a stable deep interior, since helium phase separa- 
tion from hydrogen must not necessarily occur in a metallic 
environment. Lorenzen et al. (2009, 2011) calculated the 
H-He immiscibility regions in dependence on temperature, 
pressure, and helium concentration. In metallic hydrogen, 
demixing requires a sufficiently high helium concentration 
to be energetically preferred. For instance, at 2 Mbar our 
Uranus and Neptune models have a temperature of 3800- 
4500 K. The phase diagram of Lorenzen et al. (2011) pre- 
dicts demixing under these conditions only if the He con- 
centration is above ~ 1%. Using 

T=(^-7^)", (i) 

N \ m H Y m H2 o y J 

where the m, denote the molar masses of H, He, and wa- 
ter, we calculate helium particle concentrations N-R e /N of 
1.7-3% for Z = 0.9-0.95. This is only slightly above the 
demixing condition in pure H-Hc mixtures. As the demix- 
ing behavior of helium in more complex mixtures such as 
H-HC-H2O is unknown, we cannot rule out miscibility of 
He the deep interior of Uranus and Neptune, be it stable 
(Uranus case) or unstable (Neptune case) as an explana- 
tion for the observed abundances. We encourage future 
work on the miscibility behavior of planetary mixtures, 
such as started by, e.g., Chau et al. (2011); Wilson and 
Militzer (2012a,b). 

4-3. Luminosity 

The modified shape and rotation data do not affect the 
cooling behavior. With Uranus, this is to be expected, 
since former calculations with various different structure 
models as allowed by the large error bars in J2 and J4 pro- 
duced an uncertainty of about 1 Gyr in the evolution only, 
smaller than the uncertainty induced by the observational 



error bar of T e g (Fortney et al., 2011). With Neptune, this 
is of some surprise, since the new models are ~ 1000 K 
colder in the deep interior than former models used for 
cooling curve calculations (FN10). However, the cooling 
curve depends on changes of temperature with time. These 
are similar right after the rapid cooling at young ages. 

As the heat stored in the interior after formation de- 
pends on the specific heat of the bulk material, where c v is 
smaller for rocks than for ices, future work should include 
the admixture of rocks into the envelopes and investigate 
the maximum possible shortening of Uranus' cooling time 
within the approach of homogeneous, adiabatic evolution. 

Magnetic field models of Uranus and Neptune require 
that 60-70% of the region interior to the ionic water layer, 
corresponding to 0.42-0.56 R p (Redmer et al., 2011), is 
stable to convection (Stanley and Bloxham, 2006). Using 
our interior models, this corresponds to 0.25-0.5 of the 
planet's mass. Indeed we can reproduce Uranus' measured 
luminosity if we assume that the heat flux from within 
0.45-0.5 M p is a negligible fraction of the heat flux in case 
of a convective, adiabatic interior. This consistency was 
already noted by Podolak et al. (1991). However, this 
mass or radius level is not supported by any of our Uranus 
structure models. Even if a density gradient occurs con- 
tinuously rather than as a sharp transition, our structure 
models for Uranus predict a location farther out at 0.9My. 
Also for Neptune, a stable interior up to 0.25-0.5 M p is 
not supported by our structure models, where the deepest 
strong density gradient can occur at 0.57A/n- Neptune's 
luminosity is best explained by the absence of a stable 
region. Future work should aim at finding consistent so- 
lutions for the structure, thermal evolution, and and mag- 
netic field generation of Uranus and Neptune. At the cur- 
rent stage, the envelope separation (density gradient) of 
our structure models does not seem to be directly related 
to the luminosity, nor to the magnetic field generation. 

4-4- Why a dichotomy? 

Beside the low luminosity, Uranus differs from Neptune 
in having a high obliquity, dense narrow rings and its five 
largest satellites on regular orbits, while Neptune has a 
more Saturn-like obliquity and extended dusk disk with 
diffuse rings, and two major satellites (Triton and Nereidc) 
on irregular orbits. These observed properties point to 
different formation histories and as a consequence, also 
to different internal structures. As discussed by Steven- 
son (1986), the stochastic process of impacts of various 
sizes and obliquities may have caused the differences we 
see today. A composition gradient in a heavy-element-rich 
giant planet could be a remnant of the formation process 
(Podolak et al., 1991) that has been disturbed in Nep- 
tune by a last big direct impact but not so in Uranus 
(an oblique impact). We are not the first authors to sug- 
gest a dichotomy between the internal structures of Uranus 
and Neptune. But our structure models are the first ones 
to confirm this pre- Voyager hypothesis of their different 
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structures also on the basis of computed moment of iner- 
tia values, with an absolute difference of dX ~ 0.03(~ 0.01) 
for the models with the (un-)modified shape and rotation 
data. 

4-5. A Vote for improved observational data 

New measurements of Uranus and Neptune's physical 
parameters could improve our knowledge of their bulk 
compositions and internal structures considerably. For in- 
stance, the outer envelope metallicity is rather sensitive to 
the assumed solid-body rotation period. While the modi- 
fied periods used in this paper are perhaps not the correct 
ones for Uranus and Neptune, since differential rotation is 
possible and the winds may not be in a state of minimum 
energy, the presented models clearly indicate the need for 
more data. Recently, Karkoschka (2011) analyzed a collec- 
tion of atmospheric circulation data compiled from Voy- 
ager 2 and HST observations of Neptune's atmosphere over 
a time-span of 20 years and found a rotation period close 
to the Voyager value, based on the stability of the motion 
of high-altitude clouds. Better determination of the rota- 
tion periods, shapes, Uranus' luminosity, Neptune's grav- 
ity field, and of the envelope mctallicities below the water 
cloud level are indispensable for a better understanding of 
the interior, evolution, and formation of the icy planets in 
the outer solar system. We encourage the application of 
various discovery methods in the future, and in particular, 
to fly space missions to Uranus and Neptune. 

4.6. A dichotomy in the atmosphere compositions? 

One of the indications for a dichotomy in the inte- 
rior structure is the possible difference in the outer en- 
velope mctallicities Z±, which are at most 0.08 (< 0.18) 
for Uranus but up to 0.65 (< 0.54) for Neptune, using the 
(un-)modified shape and rotation data. In case of vertical 
mixing up into the troposphere, one would expect to see 
a signature from different outer envelope metallicities in 
the lower atmosphere. Over the past 60 years, several at- 
tempts of deriving the atmosphere abundances from mea- 
sured spectra have been undertaken, see, e.g., the reviews 
of Fcgley et al. (1991); Gautier et al. (1995); Guillot and 
Gauticr (2007). However, that turned out to be quite chal- 
lenging, mostly because of the numerous absorption lines 
of CH4 and H2, uncertainties in the line shapes and ab- 
sorption coefficients, non-equilibrium chemistry, and the 
dependences of the derived abundances among each other. 
Today, carbon is the heavy element for which the interpre- 
tation of the available data gives the least ambiguous pic- 
ture. A value of C:H=30-60x solar for both planets has fi- 
nally emerged (Gautier et al., 1995). Assuming C:H= 30x 
solar and similar enrichments also of N and O, Hubbard 
et al. (1995) estimate an outer envelope ice mass fraction 
of 0.08. While this is in good agreement with the possible 
Z\ values of our Uranus and Neptune models, it does not 
support our high-Zi Neptune models as required for qual- 
itatively different interiors. On the other hand, both CO 



and HCN have been detected in the atmosphere of Nep- 
tune but not of Uranus (Gautier et al., 1995), indicating 
different processes at work in their atmospheres. To re- 
place assumptions by real data, we encourage deep entry 
probe missions to Uranus and Neptune. 

4-7. Past and present structure models 

As shown in figures 1-3, more definite conclusions about 
the dissimilarity of the internal structures are prevented by 
the current observational uncertainty in the gravitational 
moments. This is at odds with the struggles of earlier 
modelers (Podolak et al., 1995; Marlcy et al., 1995) to 
find Uranus and Neptune models at all that matched the 
less accurately known gravity data in the past. Moreover, 
despite the tighter present constraints, our models seem to 
encompass many of the previously found models. 

In particular, Hubbard and Marley (1989) developed 
Uranus models with physical EOS, smooth transitions be- 
tween an outer H-He rich envelope, an inner ice-rich enve- 
lope, and a rock core, and showed that models with 5-15% 
H-He in the inner envelope could give best agreement with 
the imposed constraints. Such an H-Hc abundance is also 
seen in our Uranus models. 

When the Voyager gravity and rotation rate data for 
Neptune became available, Podolak et al. (1995) applied 
the three-layer approach and physical EOS to compute 
Uranus and Neptune models. Their Uranus models re- 
quired some H-He to be present in the deep interior, too, 
whereas for Neptune this was found to be optional (non- 
conventional models). They could also find a different 
Neptune model with an inner envelope of pure water and 
transition deeper inside at ~ 0.7i?N, and a highly enriched 
outer envelope, indicating a possible difference in the in- 
terior structures of Uranus and Neptune. Our Neptune 
models in the upper right corners of Figs. 2,3 are similar 
to that latter one, while their non-conventional models are 
within the bulk of our Uranus and Neptune models. 

Marley et al. (1995) again allowed for smooth density 
gradients and generated models with random density dis- 
tributions. Nevertheless, they found that a rather sharp 
transition from a low-density outer envelope to a high- 
density inner envelope is necessary to fit the Voyager val- 
ues of J2 and J4, where the pressure-density relation in 
the inner envelope could be that of an ice layer. Because 
of the admittance of smooth density gradients, they found 
for the first time that the transition could occur as deep as 
between ~ 60-65% of the planet's radius. Our model N2b 
shows the same properties, although within the three-layer 
approach. 

Finally, Helled et al. (2011) were the first to use the accu- 
rate post- Voyager gravity data of Jacobson (2007, 2009). 
Their polynomial density-profiles could well represent a 
metallicity gradient that rises from about solar metallicity 
at the surface to about 85% in the center, with slightly 
higher values preferred for Neptune but no indication of 
different internal structures. The found metallicities are 
within those of our models. 
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According to the present work, the application of the 
modified shape and rotation data gives a more pronounced 
indication of different structure than seen so far. 

4-8. Conclusions 

We present the full sets of three-layer interior models 
(with H/He/water envelopes and rocks confined to the 
core) of Uranus and Neptune for different solid-body ro- 
tation periods and Battenings, using the improved gravity 
data by Jacobson (2007, 2009), and physical equations of 
state. We find that the resulting bulk composition is insen- 
sitive to the current level of uncertainty in the input data 
(observational constraints and equations of state) as our 
results arc in good agreement with previous calculations 
(e.g. Podolak et al. 1995; Marley et al. 1995; Helled et al. 
2011). However, our models with the modified rotation pe- 
riods and shapes suggest that Uranus and Neptune could 
be quite different. Uranus would have an outer envelope 
with a few times the solar metallicity which transitions to 
a heavily enriched (~ 90% by mass heavy elements) inner 
envelope at 0.9 My, giving a rather low moment of inertia 
of ~ 0.222. In Neptune, this transition can occur deeper 
inside at 0.6 Mn and be accompanied by a more moderate 
increase in metallicity, leading to a less centrally condensed 
planet with A ~ 0.255. While the observed magnetic fields 
of Uranus and Neptune are similar and can be reproduced 
by a rather narrow range of dynamo models, dissimilar in- 
teriors are required to explain the measured luminosities. 
We have presented a new indication for different internal 
structures based on the application of modified shape and 
rotation data. However, the density gradient in our mod- 
els appears to be generally farther out than required by 
evolution and magnetic field models. 

The authors thank the two referees for providing fruitful 
comments that led us to consider the results in a wider con- 
text. NN acknowledges support from the DFG RE 881/11- 
1. 
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Figure 2: (Color online) 

Heavy element mass fraction in the outer envelope (Z\) and inner envelope (Z2) of Uranus models 
(black) and Neptune models (grey) as labeled with the Voyager shape and rotation data. The 
solid lines frame the full set of solutions for each planet. Dashed lines within the box of Neptune 
models indicate solutions of same transition pressure in [Mbar] as labeled. Numbers at selected 
models (filled circles) give T c [K], P c [Mbar], M c [M®], the ice-to-tock ratio I:R, and A. The 
dotted line is a guide to the eye for the solar metallicity Zq = 0.015. 
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z 2 

Figure 3: (Color online) 

Same as Fig. 2 but using the modified shape and rotation data for Uranus (cyan) and Neptune 
(blue). Models U2, N2a, and N2b (Table 2) are highlighted by black-filled, big circles. The boxes 
of Fig. 2 are also shown to facilitate the comparison. 




Figure 4: (Color online) 



Internal density (left panel) and mass (right panel) profiles over the normalized planetary radius 
of the models Ul (solid, black), U2 (thick solid, cyan), Nl (short-dashed, grey), N2a (dotted, 
blue), and N2b (long-dashed, blue). Sec Table 2 for details. 
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Figure 5: (Color online) 

Internal density and temperature profiles over pressure. See Fig. 4 for description of line-styles 
and labels. 
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Figure 6: (Color online) 

Thermal evolution of Uranus (model U2, solid, cyan) assuming that the heat flux from the mass 
interior to the cut-off mass (x-axis) is limited to a fraction / of the convective energy flux. For 
/ = 0-0.4 and a cut-off mass of 0.45-0.85, a cooling time in agreement with the age of the 
solar system (horizontal dashed line) can be found. For comparison, we also show the same 
calculations for Neptune (model N2a, dotted, blue). 
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Table 1: Physical constraints of Uranus and Neptune. 



Parameter 


U data 1 


U data 2 


N data 1 


N data 2 


M p (M e ) 


14.536 


14.536 


17.148 


17.148 


w/2?r 


17h 14m 40s a 


16h 34m 24s c 


16h 6m 40s a 


17h 27m 29s c 


Ti (K) 
Toff (K) 


76(2) b 
59.1(3) s 


76(2) 
59.1(3) 


72(2) a 
59.3(8) s 


72(2) 
59.3(8) 


^cq 

J2/10" 2 
4/10- 4 


26,200 d ' c 


26,200 d <° 


25,225 d -' 


25,225 d <' 


0.334129(72) d < c 
-0.3044(102) d 


0.334129(72) d ' c 
-0.3044(102) 


0.340843(450) d ' f 
-0.334(29) d - f 


0.340843(450) d ' f 
-0.334(29) d ' f 


R cq (km) 

J2/IO- 2 

J4/IO- 4 


25,559(4) a 
0.351099(72) c 
-0.3361(100) c 


25,559(4) c 
0.351099(72) c 
-0.3361(100) c 


24,766(15) a 
0.35294(45) c 
-0.358(29) c 


24,787(4) c 
0.35294(45) c 
-0.358(29) c 



Numbers in parenthesis are the observational error bars in the last digits. The gravitational moments J' 2n arc the 

measured ones. They refer to a reference equatorial radius R' cq . The gravitational moments Jin refer to the equatorial 

radius at the 1-bar pressure level, i? e q- 

a Lindal (1992), 

b Lindal ct al. (1987), 

c Helled et al. (2010), 

d http: / /ssd.jpl. nasa.gov, 

°Brozovic and Jacobson (2009), 

f Jacobson (2009) 

b after Guillot and Gautier (2007) 
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Table 2: Five resulting structure models. 



Parameter 


Ul 


U2 


Nl 


N2a 


N2b 


R cq (km) 


25,559 


25,559 


24,773 


24,786 


24,786 


-Rmean (km) 


25,388 


25,378 


24,622 


24,650 


24,650 


/ 


0.0198 


0.0210 


0.0180 


0.0163 


0.0163 


J2/10- 2 


0.35107 


0.35107 


0.3533 


0.3531 


0.3530 


J4/10- 4 


-0.345 


-0.344 


-0.378 


-0.383 


0.3835 


Pi _ 2 (Mbar) 


0.15 


0.15 


0.10 


0.30 


1.50 


mi_ 2 (M p ) 


0.913 


0.918 


0.953 


0.864 


0.571 


ri-2 (-Rp) 


0.772 


0.757 


0.927 


0.823 


0.605 


Zi 


0.17 


0.08 


0.30 


0.60 


0.644 


z 2 


0.915 


0.944 


0.833 


0.805 


1 


M c (M e ) 


0.61 


0.36 


3.15 


3.02 


0.35 


To (K) 


6000 


6500 


5500 


5000 


5200 


P c (Mbar) 


5.5 


6.0 


6.0 


5.5 


8.3 


P (Mbar) 


8.3 


7.9 


16.4 


15.3 


10.5 


M z (M e ) 


12.4 


12.7 


14.4 


13.9 


14.4 


I:R 


19.2 


35.3 


3.7 


3.6 


13.7 


r (Gyr) 


10.0 


9.1 


4.8 


4.4 


4.3 


A 


0.2296 


0.2224 


0.2405 


0.2555 


0.2557 



Uranus model Ul and Neptune model Nl are based on the input data sets No. 1, while models U2, N2a, and N2b are 
on data sets No. 2, see Table 1. The pressure Pn refers to m = 0, whereas P c to m = M c . 
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Table 3: Example of a tabulated interior profile. 



Mass 


Pressure 


Radius 


Temperature 


Density 




(GPa) 




(Rm) 


(K) 


(g/ccm) 


14.5322753 


1.0000E-04 


3 


.979310 


76.0 


4.4876E-04 


13.2697748 


1.4999E+01 


3 


.073483 


2338.8 


4.0541E-01 


13.2697239 


1.5000E+01 


3 


.073465 


2338.8 


1.1879E+00 


0.5987030 


5.5226E+02 





.700230 


6083.5 


4.0700E+00 


0.5984897 


5.5234E+02 





.700034 


6083.6 


9.0790E+00 


0.0000019 


8.2114E+02 





.010000 


6083.6 


1.0326E+01 



Tables of the interior profiles of the models Ul, U2, Nl, and N2b arc published in their entirety as supplemental material 
in the electronic edition. A portion of the table of model Ul is shown here for guidance of the tables' form and content. 
The selected rows show, from top to bottom, the outer boundary, the outer/inner envelope transition, the core- mantle 
boundary, and the center. 
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